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Abstract 

We use hybrid functionals and restricted self-consistent GW, state-of-the-art theoretical ap- 
proaches for quasiparticle band structures, to study the electronic states of delafossite Cu(Al,In)02, 
the first p-type and bipolar transparent conductive oxides. We show that self-consistent GW gives 
remarkably wider band gaps than all the other approaches used so far. Accounting for polaronic 
effects in the GW scheme we recover a very nice agreement with experiments. Furthermore, the 
modifications with respect to the Kohn-Sham bands are strongly fc-dependent, which makes ques- 
tionable the common practice of using a scissor operator. Finally, our results support the view 
that the low energy structures found in optical experiments, and initially attributed to an indirect 
transition, are due to intrinsic defects in the samples. 

PACS numbers: 71.20.-b 71.45. Gm 78.20.-e 71.15.Qe 71.35. Cc 
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Many high-technology devices, such as flat panel displays, touch screens, or even thin- 
film solar cells, require the use of thin transparent contacts. These contacts are usually 
built from insulating oxides that, for a certain range of doping, become conductive while 
retaining transparency in the visible spectrum. The most common examples of these so- 
called transparent conductive oxides (TCOs) are electron (n-)doped SnC>2, 1^03, and ZnO. 
Hole (p-)doping of wide gap semiconductors was for long time very hard to obtain [l]. It is 
therefore not surprising that the discovery of p-doping in CuA102 thin films with a carrier 
mobility of about 10cm 2 /(Vs) attracted great interest [2J. Other members of the delafossite 
family, like CuGa02 [3] and Culn02 Q], were discovered shortly after. The latter compound 
is particularly interesting as it exhibits bipolar (n- and p-type) conductivity by doping with 
appropriate impurities and tuning the film- deposit ion conditions [4j. This opens the way to 
the development of transparent p-n junctions, and therefore fully transparent optoelectronic 
devices, functional windows and stacked solar cells with improved efficiency. 

CuA102 is by far the most studied system of the family of delafossite TCOs. However, 
there is still no agreement either on the origin of the p-type conductivity, or on the electronic 
bands of the pure crystal. Measurements of the direct optical band gap (E$ ir ) of CUAIO2 

h fin g 

fall in the range from 2.9 to 3.9 eV |2|, |5H8], with most values in the interval 3.4-3.7eV. 
These experiments also yield a large dispersion of indirect gaps (E 1 ^), from 1.65 to 2.1 eV, 



with one experiment measuring 2.99 eV [8J. Unfortunately, there is only one photoemission 
experiment j^] that yields 3.5 eV for the quasiparticle band gap. Note that the optical 
and quasiparticle gaps differ by the exciton binding ener gy. Concerning Culn02, optical 



experiments measured E^ between 3.9 and 4.45 eV 



101 ] . with only one estimation of 



io|. 



Eg* at 1.44 eV 

From the theoretical perspective, the situation is also quite complex, even if the full Cu 
3d shell should exclude the strongly correlated electron regime. These materials are usually 
studied within density functional theory (DFT), using the standard local density (LDA) or 
generalized gradient approximations (GGA). However, it is well known that the Kohn-Sham 
band structures systematically underestimate the band gaps. For similar compounds, like 
CU2O and CuIn(S,Se)2, Kohn-Sham LDA calculations lead to unreasonable band structures, 
in particular due to the misrepresentation of the hybridization between the d electrons of the 
metal and p electrons of the anion 

HQ . To overcome this situation, hybrid functionals 



have been recently proposed, with very promising results L3[, especially for materials with 
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small and intermediate band gaps 23, |l5j. Other approaches include LDA+U, that tries to 



improve the description of correlation through the introduction of a mean-field Hubbard-like 
term. This method has been quite successful in the study of strongly correlated systems, 
but it relies on a parameter U, that is often adjusted to experiments. 

Arguably the most reliable and used ab initio technique to obtain quasi-particle band 
structures is the many-body GW approach [16|]. The common practice within this framework 
is to start from a DFT calculation, and evaluate perturbatively the GW energy corrections 
to the band structure. This procedure, which we will refer to as G W , is justified when 
the departure wave functions and band structure are already close to the quasiparticle 
ones. This is indeed the case in many systems, explaining why Go Wo has been extremely 
successful in describing electron addition and removal energies for metals, semiconductors 
and insulator s | l7| . However, it has been recently shown that Go Wo fails for many transition 



metal oxides 111 



To solve this problem one can perform restricted self-consistent (sc) GW 19] . This 
technique has the advantage of being independent of the starting point at the price of 
large computational complexity. Fortunately, there is an alternative procedure that yields 



wavefunctions that are extremely c 
sc-COHSEX as explained in Ref. 



ose to those obtained in a full sc-GW calculation, namely 



20] . The dynamical effects that are absent in COHSEX 



calculations can then be accounted for by performing a final perturbative GW step. This 



method, that we will refer to as sc-GW, has been applied to many oxide compounds, 
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20]. 
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excellent results for the band gaps and the quasiparticle band structure 

Note that these theoretical techniques yield quasiparticle bands, and not optical gaps. 
To evaluate these latter quantities one mostly resorts to the solution of the Bethe-Salpeter 
equation. For the delafossite structures there is one such calculation starting from a GGA+U 
band structure [211 ] . It yields for CuA10 2 a very large exciton binding energy of about 0.5 eV 



for the first direct transition. The choice of U was found to have strong consequences on 
the width of the band gap, but it did not affect significantly the exciton binding energy. We 
can thus assume that 0.5 eV is a reasonable estimate. 

In the following, we present calculations of the band structures of CuA102 and Culn02 
using some of the most accurate theoretical tools available in the community. These include 
the standard LDA, hybrid functionals (namely B3LYP 22j and two flavors of Heyd-Scuseria- 
Ernzerhof, HSE03 and HSE06 [23]), LDA+U, G W and sc-GW. As discussed above, we 
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FIG. 1: (Color online) Band gaps of CuA10 2 using: LDA, LDA+.7, hybrid, G Q W Q , sc-GW, and 
sc-GW including model polaronic corrections. The horizontal zones contain data extracted from 
various optical experiments (see text). 

expect sc-GW to be the most accurate ab initio approach. When the comparison was 
possible, we found our results in excellent agreement with previous calculations (Refs. (5, 7, 
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24 



26| for LDA, Ref. |26| for B3LYP, and Ref. |27j for GGA+U). 
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and ABINIT [29 1 



The hybrid and LDA+U calculations were performed with VASP 
respectively, using the PAW formalism and an energy cutoff of 44 Ha. The parameter U was 
set to 8eV as in Ref. 21] . Our GW calculations were performed with ABINIT, starting from 



LDA band structures and using norm-conserving pseudopotentials with semicore states (3s 
and 3p for Cu and 4s and 4p for In) included in the valence. The energy cutoff was 120 Ha 
for the ground state calculation, and the k-point grid was a 4 x 4 x 4 Monkhorst-Pack. As 
the experimental and LDA relaxed geometries are very close (within 1%), and the small 
contraction of the lattice in LDA has a ne glig ible effect on band structures (< 0.05 eV), we 



employed experimental lattice parameters 



261 ] . Note that it was absolutely essential to use 
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FIG. 2: (Color online) Band structures for CuAlC^: comparison of LDA (red dashed lines) with 
sc-GW (left panel), HSE03 (central panel), and LDA+U (right panel). 



301 ] . due to the extremely slow convergence with respect to the number 



the method of Ref. 
of empty states. 

In Figs. [T]and [2] we show direct and indirect photoemission gaps and the band structures 
of CuAlC>2 obtained using different theoretical approaches. The minimum E du of CuAlOo is 



always found at L, where the dipole transition between the band edge states is allowed [25 ]. 
All calculations, except sc-GW, give a fundamental E^ d between the conduction band min- 
imum at T and the valence band maximum along the T-F line. The experimental data for 
optical gaps are also presented with an error bar that reflects the dispersion of the most 
likely values found in literature. LDA exhibits, as expected, the smallest gaps. Basically 
every approach beyond it opens up the gap by different amounts and modifies the band 
dispersions. The direct and indirect gaps have similar behaviors in the different theories, 
and both increase when going from LDA<G W <HSE03<HSE06<B3LYP<sc-G'W. On 
the other hand, the difference E dn — E^ d seems to decrease with the sophistication of the 
method, reaching nearly zero for the sc-GW calculation. This is a consequence of the drastic 
change of the conduction band dispersion, which displaces the conduction minimum from T 
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FIG. 3: (Color online) Band gaps of Culn0 2 using: LDA, LDA+J7, hybrid, G W , sc-GW. The 
horizontal zones contain data extracted from optical experiments (see text). 

to L when sc-GW is applied (see Fig. |2]). Only LDA+[7 does not follow the trend, as it is 
the only case in which E dir — E™ d gets significantly larger than in LDA. 

Looking at the direct gap, we point out that most of the methods give results that are 



within the experimental range, when an exciton binding energy of around 0.5 eV [2l| is 
considered. This is true for LDA+U, GqWq, the hybrids HSE03 and HSE06. However, 
for sc-GW and even for B3LYP, the theoretical gap is larger by about 1-1.5 eV than the 
experimental findings. For Culn02 (see Fig. [3]) we have to make the comparison with care, as 



. A meaningful 



the smallest E dn is located at T, where optical transitions are forbidden 
comparison with experiments must consider the gap at L. Thus, we find that both trends 
and quantitative results are analogous to those for CuA10 2 . In particular, sc-GW yields 
again E dir larger by 1-1.5 eV than the experimental range. 

We stress again that, to date, sc-GW is arguably the best method available to estimate 
band gaps of wide-gap semiconductors, and that it gives excellent results for compounds like 



CU2O and CuIn(S,Se)2 [III, ll2j. It is unlikely that the presence of defects can lead to such 
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a large shrinkage of E dn . However, there is another effect that has been neglected up to 



now: the change of screening due to the polarization of the lattice 



31] . In fact, according to 



the experimental data [32], unfortunately available only for CuAlC>2, the polaron constant 
for this system is large (a p ~ 1), indicating a non-negligible contribution of the lattice 
polarization to the electronic screening. In other ionic compounds with similar polaron 
constants this can lead to a shrinkage to the band gap by about leV 33]. A full sc-GW 
calculation including in an ah initio framework the effects of the lattice polarization is to 
date beyond reach. However, a reliable estimate can be obtained using the model proposed 
by Bechstedt et al. [33(, which gives a static representation of the polaronic effects based 
on difference of experimental static dielectric constants. By performing a perturbative GW 
step including model polaronic effects on top of the sc-COHSEX, we found a uniform (It- 
independent) shrinkage of the band gap by 1.2 eV. As we can see in Fig. [TJ this correction 
brings our results for E dlT well within the experimental range (once the excitonic correction 



of about 0.5 eV is also considered). As it is observed in Ref. [33(, the polaronic model 
employed can only overestimate the correction. All these results point to the conclusion 
that the agreement of the other methods with experiment was fortuitous and due to a 
cancellation of errors. 

Looking now at the indirect gap, we focus on Fig. [T] as there are more experimental data 
for CuA102- All the hybrids and GW calculations yield indirect gaps much larger than 
the experimental range 1.65-2.1 eV, even taking into account any possible excitonic and 
polaronic effects. Moreover, sc-GW, the best method used in this work, yields the highest 
E™ d at around 5eV, while the difference Eg* - E™ d is in general much smaller than the 
experimental value (~2eV), and even vanishing for the sc-GW calculation. From Fig. [3] we 
realize that these conclusions are as well valid for Culn02, where the best estimates for the 
indirect band gap is much larger than the experimental value of 1.44 eV lo| . 

These are very strong arguments in favor of Robertson et al. 26] that suggested that 
the experimental "indirect gap" absorption was due to defects, and should not be present in 
the defect-free compound. Also Pellicer- Porres et al. js| questioned the interpretation of 
the low energy peaks as indirect transitions, as the absorption coefficient is more than two 
orders of magnitude larger than in typical indirect absorption edges. The most promising 
defects are oxygen interstitials O;, as LDA calculations predict low formation energies and 
the introduction of states in the gap at 0.7 and 1.4 eV [34]. However, a full clarification 
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of this issue will require sc-GW or hybrid calculations for these, and other more complex 
defects. 

Finally, we analyze more in detail the band structures of CuA102 shown in Fig. [2j LDA 
calculations (red dashed lines) are compared with sc-GW, HSE03, and LDA+t/ calculations. 
The main effect of LDA+[7 is to open the LDA gap by an amount that can be controlled 
by the parameter U. The difference E da — E^ d is in this approximation enhanced, due to 
a change of the character of the lowest conduction band along the symmetry lines. Hybrid 
calculations using HSE03 give a comparable E dir and a modified dispersion of both valence 
and conduction states close to the Fermi energy, which reduces E dir — E^ d . The conduction 
band minimum (CBM) within HSE03 is still located at T, but the difference between the 
CBM at L and T gets significantly smaller. For sc-GW, besides the further increase of the 
band gaps, the dispersion of the bands is strongly affected by the many-body effects. In 
fact, the GW corrections exhibit an unusual dispersion of around leV when looking at the 
different k-points, displacing the CBM from T to L. We note that often in semiconductor 
physics one assumes that the quasiparticle corrections can be modeled by a rigid shift (the 
so-called scissor operator). From our results it follows that one should refrain from using 
this simple approximation for these important materials. We can also conclude that hybrid 
calculations give a better description of band dispersions than LDA+J7, even if the two 
approaches yield similar band gaps. 

In conclusion, it is clear that the delafossite family exhibits complex and unusual band 
gap physics that can not be captured by standard theoretical approximations. We found that 
the direct band gap is well reproduced by the best many-body approaches if polaronic effects 
are taken into account. We can expect that this situation, of a large gap that is reduced 
substantially by polaronic effects, is quite general and is present in many more materials 
that previously expected. In fact, the apparent good agreement between calculated gaps 
(with hybrid functionals or GqWq) and experimental gaps for materials as simple and widely 
studied as LiF can be accidental, as preliminary calculations confirm: the underestimation 
of the gap by these methods (the scGW gap is indeed 2eV larger than the experimental 
and GqWo gap) is compensated by the neglect of large polaronic effects. Furthermore, the 
modifications with respect to the LDA Kohn-Sham bands are strongly fc-dependent, which 
makes questionable the common practice of using a scissor operator. The band dispersion 
obtained by hybrid functional calculations is in between the LDA and sc-GWdispersion, 
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while the LDA+U calculations open up the gap but do not give a significant improvement of 
the band dispersion. Finally, our calculations rule out the interpretation of the low energy 
features in the absorption spectra as arising from a putative indirect band gap. These 
structures should rather come from intrinsic defects, as proposed in Refs. 26]. However, 
a complete understanding of the electronic and excitation properties of these systems will 
only be achieved, in our opinion, by a high-level theoretical scheme (like sc-GW) including 
defects and effects from the lattice polarization in an ab initio framework. Work along these 
lines is in progress. 
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